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binary-source interpretations. 


the 3L1S analysis of the KMT-2021-BLG-1122. 


1. Introduction 


Light curves of gravitational microlensing events often exhibit 
deviations from the smooth and symmetric form of the event 
involved with a single lens mass and a single source star 
(Paczyński 1986). The most common cause of such deviations 
is the binarity of the lens, 2L1S event (Mao & Paczynski 1991), 
or the source, 1L2S event (Griest & Hu 1993; Han & Gould 
1997). Anomalous lensing events comprise about 10% of 
more than 3000 lensing events that are being annually de- 
tected from the three currently working lensing surveys con- 
ducted by the Optical Gravitational Lensing Experiment (OGLE: 
Udalski et al. 2015), Microlensing Observations in Astrophysics 
(MOA: Bond et al. 2001), and Korea Microlensing Telescope 
Network (KMTNet: Kim et al. 2016) groups. The prime goal 
of these survey experiments is finding extrasolar planets, which 
are detected via anomalous signals in lensing light curves. Plan- 
etary microlensing signals appear in various shapes and often 
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ABSTRACT 


Aims. We systematically inspect the microlensing data acquired by the KMTNet survey during the previous seasons in order to 
find anomalous lensing events for which the anomalies in the lensing light curves cannot be explained by the usual binary-lens or 


Methods. From the inspection, we find that interpreting the three lensing events OGLE-2018-BLG-0584, KMT-2018-BLG-2119, and 
KMT-2021-BLG-1122 requires four-body (lens+source) models, in which either both the lens and source are binaries (2L2S event) 
or the lens is a triple system (3L1S event). Following the analyses of the 2L2S events presented in Han et al. (2023), here we present 


Results. It is found that the lens of the event KMT-2021-BLG-1122 is composed of three masses, in which the projected separations 
(normalized to the angular Einstein radius) and mass ratios between the lens companions and the primary are (52, q2) ~ (1.4,0.53) 
and (53,q3) ~ (1.6,0.24). By conducting a Bayesian analysis, we estimate that the masses of the individual lens components 
are (Mı, M2, M3) ~ (0.47 Mo, 0.24 Mo, 0.11 Mo). The companions are separated in projection from the primary by (a12,a13) ~ 
(3.5, 4.0) AU. The lens of KMT-2018-BLG-2119 is the first triple stellar system detected via microlensing. 


Key words. Gravitational lensing: micro — (Stars:) binaries: general 


confused with anomalies produced by other causes. In order to 
sort out planetary signals from those of other origins, therefore, 
it is essential to analyze all anomalous events. 

In the current microlensing experiments, anomalous events 
are being analyzed almost in real time with the progress 
of events. This became possible with the implementation 
of early warning systems: the OGLE early warning sys- 
tem (Udalski etal. 1994), MOA real-time analysis system 
(Bond et al. 2001), and KMTNet EventFinder AlertFinder 
(Kim et al. 2018b) system.! These systems enable one to de- 
tect lensing events in their early phases and to promptly iden- 
tify anomalies appearing in lensing light curves. Analyses of 
the events with identified anomalies are conducted by multi- 


' Tn addition, the KMTNet EventFinder system (Kim et al. 2018a) car- 
ries out end-of-season analyses that recover many events that were 
missed for one reason or another by AlertFinder (Kim et al. 2018b) sys- 
tem. 
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ple modelers of the individual survey groups as the anomalies 
proceed, and models of the events found from these analyses 
are circulated or posted on web pages? to inform researchers in 
the microlensing community of the nature of the anomalies. In 
most cases, anomalous events are interpreted with a binary-lens 
(2L1S) or a binary-source (1L2S) model. 

For a minor fraction of anomalous events, it is known that 
the anomalies cannot be explained by the usual 2L1S or 1L2S 
interpretation. Han et al. (2023), hereafter Paper I, conducted a 
systematic investigation of events found from the KMTNet sur- 
vey during the previous seasons in search for anomalous events 
for which no plausible models had been presented. From this in- 
vestigation, they found that the anomalies appearing in the light 
curves of the events OGLE-2018-BLG-0584 and KMT-2018- 
BLG-2119 required a four-body (lens+source) model, in which 
both the lens and source are binaries (2L2S event). In this pa- 
per, we present the analysis of another four-body lensing event 
KMT-2021-BLG-1122. Unlike the two 2L2S events presented in 
Paper I, the lens system of KMT-2021-BLG-1122 is composed 
of three masses (3L1S). 

We present the analysis of the event KMT-2021-BLG-1122 
according to the following organization. In Sect. 2, we depict 
the observations conducted to obtain the photometric data of the 
event, and describe the characteristic features of the anomaly 
appearing in the light curve of the event. In Sect. 3, we begin 
by describing the lensing parameters used in the modeling con- 
ducted under various configurations of the lens and source sys- 
tem, and we then describe detail of the analyses conducted un- 
der the individual configurations in the following subsections: 
2L1S model in Sect. 3.1, 2L2S model in Sect. 3.2, and 3L1S 
model in Sect. 3.3. In Sect. 4, we specify the source star of the 
event and estimate the angular Einstein radius of the lens sys- 
tem. In Sect. 5, we estimate the physical parameters including 
the masses of the individual lens components, the distance to 
the lens system, and the projected separations among the lens 
components. We summarize the results from the analysis and 
conclude in Sect. 6. 


2. Lensing light curve and observations 


The lensing event KMT-2021-BLG-1122 was detected solely by 
the KMTNet survey. The event occurred on a faint source star, 
with a baseline magnitude of [pase = 20.36, located toward the 
Galactic bulge field with equatorial coordinates (RA, DEC) = 
(17:35:51.10, —28:26:47.22), which correspond to the Galac- 
tic coordinates (l, b) = (—0°720, 2°076). The source lies in the 
KMTNet BLG14 field, toward which observations were con- 
ducted with a 1 hr cadence. The event was identified by the 
KMTNet AlertFinder system and an alert was issued on 2021 
June 3, which corresponds to the abridged Heliocentric Julian 
Date of HJD’ = HJD — 2450000 ~ 9368.5. 

Observations of the event were carried out with the use of 
the three 1.6 m telescopes operated by the KMTNet group. The 
telescopes are distributed in the three continents of the Southern 
Hemisphere for continuous coverage of lensing events, and the 
sites of the individual telescopes are the Siding Spring Observa- 
tory in Australia (KMTA), the Cerro Tololo Interamerican Ob- 
servatory in Chile (KMTC), and the South African Astronomical 
Observatory in South Africa (KMTS). The camera mounted on 
each telescope is composed of four 9k x 9k CCD chips, which in 
combination yield a 4 deg? field of view. 


2 For example, the model page of lensing events maintained by 
Cheongho Han (http: //astroph. chungbuk.ac.kr/~cheongho). 
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Fig. 1. Light curve of KMT-2021-BLG-1122 constructed with the 
combined data obtained from the three KMTNet telescopes (KMTA, 
KMTC, and KMTS). The curve drawn over the data points is the best- 
fit model found under the 2L1S interpretation. The lower panel displays 
the whole view, and the upper panel shows the enlarged view of the 
anomaly region. The epochs marked by tı, t2, and t, indicate the three 
major anomaly features. 


Images containing the source were mainly taken in the 7 
band, and about one tenth of images were obtained in the V band. 
Reductions of the images and photometry of the event were done 
using the KMTNet pipeline constructed based on the pySIS code 
of Albrow et al. (2009), and additional photometry was carried 
out for a subset of the KMTC data set using the pyDIA code of 
Albrow (2017). The analysis of the event was done based on the 
I-band light curve constructed from the pySIS reduction, and the 
pyDIA photometry data were used for the source color measure- 
ment. We discuss the detailed procedure of the source color mea- 
surement in Sect. 4. According to recipe described in Yee et al. 
(2012), we readjusted the error bars of the data determined by 
the automated pipeline in order that the errors were consistent 
with the scatter of data and y? per degree of freedom (dof) for 
each data set became unity. 


Figure | shows the lensing light curve of KMT-2021-BLG- 
1122 constructed with the combined /-band pySIS photome- 
try data obtained from the three KMTNet telescopes. It shows 
a strong anomaly relative to a 1L1S form, and the anomaly 
displays a complex pattern that is characterized by the three 
strong features appearing at HJD’ ~ 9369.7 (ti), ~ 9371.7 (t2), 
and ~ 9373.3 (t3). From the characteristic profiles described in 
Schneider & Weiss (1986), the anomalies at around f, and f3 ap- 
pear to be a pair of caustic-crossing features produced when 
a source enters and exits a caustic, respectively. On the other 
hand, the profile at around t2, which is approximately symmet- 
ric with respect to the peak, appears to be produced when a 
source approaches close to or crosses over the cusp of a caus- 
tic (Schneider & Weiss 1987). 
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Fig. 2. Enlarged view around the central anomaly region of the light curve. Three model curves obtained under the 3L1S, 2L2S, and 2L1S lens- 
system configurations are drawn over the data points, and the three lower panels show the residuals from the individual models. The two insets in 
the top panel show the lens-system configurations of the 2L1S and 2L2S models. In each inset, the red figure is the caustic and the line with an 
arrow represents the source trajectory. For the 2L2S configuration, there are two source trajectories, in which the black and blue ones represent the 


trajectories of the primary and secondary source stars, respectively. 


In general, the light curve profile of a 2L1S event in the 
region between a pair caustic-crossing spikes exhibits a "U"- 
shape pattern. Occasionally, deviations from this U-shape pat- 
tern can arise if a source sweeps one fold of a caustic, for ex- 
ample, OGLE-2016-BLG-0890 (Han et al. 2022c), but, in this 
case, the resulting deviation is smooth and thus differs from the 
sharp pattern of the observed anomaly at around t2. This sug- 
gests that a different explanation is needed to explain the origin 
of the anomaly. 


3. Interpreting the anomaly 


In order to explain the anomalies in the lensing light curve, we 
conducted a series of modeling under various interpretations of 
the lens-system configuration. We tested three models under the 
2L1S, 2L2S, and 3L1S configurations. We exclude 1L2S and 
1L3S models because the observed light curve exhibits obvious 
caustic-crossing features, that is, those around t, and t3, while 
1L2S and 1L3S events do not involve caustics.’ 


3 There has been only one 1L3S microlensing event published to date: 
OGLE-2015-BLG-1459 (Hwang et al. 2018). 


In the modeling under each tested interpretation, we search 
for a lensing solution, representing a set of lensing parameters 
that characterize the observed lensing light curve. The basic lens- 
ing parameters of a 1L1S event include (tọ, uo, tg), which denote 
the time of the closest source approach to the lens, the source- 
lens separation (impact parameter) at fo, and the Einstein time 
scale, respectively. The Einstein time scale represents the time 
required for the source to transit the angular Einstein radius of 
the lens, 6g, and the length of uo is scaled to 6g. 


A 2LI1S system corresponds to the case in which the lens 
contains an extra component compared to the 1L1S system. Con- 
sideration of the extra lens component requires one to include 
additional parameters in modeling. These additional parameters 
are (s, q, œ, p), and they denote the projected separation (normal- 
ized to 6g), and mass ratio between the lens components (Mı 
and M2), the angle (source trajectory angle) between the relative 
lens-source proper motion, jl, and the axis connecting Mı and 
Mp, and the ratio between the angular radius of the source, 0+, 
to the Einstein radius, that is, p = 0,/@¢ (normalized source ra- 
dius), respectively. The normalized source radius is included in 
modeling because a 2L1S event usually involves caustic cross- 
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ings, during which lensing magnifications are affected by finite- 
source effects (Bennett & Rhie 1996). For the computation of 
finite-source magnifications, we utilize the map-making method 
of Dong et al. (2008). 


For the modeling of a 2L2S system, it is required to in- 
clude extra parameters in addition to those of the 2L1S model 
in order to describe an extra source. These extra parameters in- 
clude (to.2, u0,2, P2, qF), Where the first three denote the closest 
approach time, impact parameter, and normalized radius of the 
source companion, and the last parameter indicates the flux ratio 
between the companion (S2) and primary (S1) of the source. In 
the 2L2S modeling, we use the notations (fo,1, Uo,1, p1) to denote 
the parameters related to S1. 


A 3LIS system has an extra lens component (M3) compared 
to the 2L1S system, and this requires one to include additional 
parameters in modeling. These parameters include (53, q3, Y), 
which denote the separation and mass ratio between Mı and M3, 
and the orientation angle of M3 as measured from the M,—M) 
axis with a center at the position of M4, respectively. In order to 
distinguish the parameters describing M> from those describing 
M3, we use the notations (s2,g2) to denote the separation and 
mass ratio of M2. The lensing parameters for the 1L1S, 2L1S, 
2L2S, and 3L1S models are summarized in Table 1 of Paper I. 


3.1. 2L1S interpretation 


We started with modeling the light curve under the 2L1S inter- 
pretation of the anomalies. In this modeling, we searched for the 
binary parameters s and q via a grid approach, while we seek 
the other parameters via a downhill approach using a Markov 
Chain Monte Carlo (MCMC) algorithm. We first constructed a 
Ay? map on the plane of the grid parameters, identified local 
minima, and then refined the individual minima by allowing all 
lensing parameters to vary. From the 2L1S modeling, we found 
no lensing solution that could describe all the anomaly features, 
and this explains why the event had been left without a plausible 
model. 


Although the light curve could not be explained by a 2L1S 
model, we found a solution that could approximately describe 
the caustic-crossing features at around t, and t3. The model curve 
of this solution, found from modeling the light curve with the 
exclusion of the data around t, is drawn over the data points 
in Figures | and 2, and the corresponding configuration of the 
lens system is presented in the upper inset of the top panel 
of Figure 2. The binary lensing parameters of the solution are 
(s,q) ~ (1.69, 0.43), which results in a single resonant caustic 
elongated along the binary axis. According to this solution, the 
source entered the caustic by crossing the caustic fold lying just 
above the left side on-axis cusp, passed through the inner region 
of the caustic, and then exited the caustic by crossing the upper 
fold of the caustic. The caustic entrance and exit resulted in the 
anomaly features matching well the observed ones around ft, and 
t3, respectively. 


3.2. 2L2S interpretation 


The fact that the two caustic-crossing features can be described 
by a 2L1S model suggests that the other anomaly feature around 
ty may be explained with the introduction of an extra source 
(2L2S) or an extra lens component (3L1S) to the 2L1S configu- 
ration. In this section, we present the analysis conducted under 
the 2L2S interpretation. 
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Table 1. Lensing parameters of the 2L2S solutions 


Parameter Static Orbital 
x /dof 1729.1/1426 1453.6/1429 
to, (HJD’) 9371.208 + 0.026 9370.681 + 0.007 
uo2 0.057 + 0.002 0.044 + 0.001 
to2 (HJD’) 9370.810 + 0.024 - 
Uo2 0.194 + 0.004 — 
tg (days) 14.36 + 0.31 14.87 + 0.11 
sS 1.723 + 0.020 1.447 + 0.006 
q 0.649 + 0.043 0.261 + 0.002 
æ (rad) 2.615 + 0.009 2.415 + 0.005 
Ss - 0.079 + 0.001 
qs = 0.103 + 0.001 
ġ (rad) - 1.326 + 0.012 
ds,/dt (yr) - —0.064 + 0.031 
d¢/dt (deg/day) - 53.79 + 0.30 
pı (1073) 3.29 + 0.17 4.58 + 0.05 
p (1077) 2.23 + 0.79 0.10 + 0.01 
qF 0.71 + 0.05 0.15 + 0.01 


Notes. HJD’ = HJD — 2450000. 


The 2L2S modeling was carried out based on the 2L1S solu- 
tion. With the initial parameters obtained from the 2L1S model- 
ing, we searched for a 2L2S solution by testing various trajecto- 
ries of the second source considering the location and magnitude 
of the anomaly around fj. We first check the case in which the 
position of the source companion with respect to the primary 
source does not vary: "static 2L2S model". The model curve of 
the best-fit static 2L2S solution is presented in Figure 2, and the 
lensing parameters are listed in Table 1. The lens-system con- 
figuration of the solution is shown in the middle inset of the top 
panel. The caustic configuration is similar to that of the 2L1S so- 
lution except that there are two source trajectories of Sı and S2, 
which are marked by black and blue arrowed lines, respectively. 
According to the solution, the source companion approached 
the lens slightly earlier than the primary source with an impact 
parameter greater than that of the primary source. The second 
source passed the tip of the upper left caustic cusp, and this pro- 
duced the anomaly feature around t2. Although the static 2L2S 
model can give rise to a sharp anomaly features appearing in the 
inner region between the caustic spikes, it was found that the 
residual from the 2L2S model was substantial, as shown in the 
second lower panel of Figure 2. 


We further check whether the anomalies can be explained 
by considering the orbital motion of the source: "orbital 2L2S 
model". We consider the source orbital motion by introducing 
five extra parameters of (Ss, gs, 6, ds;/dt, db/dt), where sẹ, and 
qs represent the normalized separations and mass ratio between 
the source components, respectively, ¢ is the orientation angle of 
S2 with respect to S1, and (ds;/dt, dé/dt) represent the change 
rates of the s, and ¢ induced by the source orbital motion, re- 
spectively. It is found that the consideration of the source orbital 
motion substantially improves the fit, by Ay? = 275.5, with re- 
spect to the static model. The model curve of the best-fit orbital 
2L2S solution is drawn in the top panel of Figure 2, residual 
is shown in the panel labeled "2L2S (orbit)", and the lensing 
parameters are listed in Table 1. According to the model, the 
anomaly features around ¢; and t3 were produced by the pri- 
mary source as is in the static case, and the feature around fz was 
produced by the secondary source. We note that the secondary 
source crossed the caustic two more times at HJD’ ~ 9368.7 (to) 
and ~ 9372.7, and the former crossing is not covered by the data 
and the latter crossing feature is difficult to identify in the data. 
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Fig. 3. Map of Ay” on the log s3-log q3 parameter plane obtained 
from the initial grid searches for (s3, q3) parameters. Color coding is set 
to represent points with < Ino (red), < 2no (yellow), < 3no (green), 
< 4no (cyan), < 5no (blue), and < 6no (purple), where n = 9. The 
upper panel shows the map around the minima constructed with denser 
grids. 


Table 2. Lensing parameters of the best-fit 3L1S solution 


Parameter Value 
x /dof 1399.8/1425 
to (HJD’) 9370.226 + 0.061 
Ug 0.126 + 0.008 
te (days) 14.74 + 0.93 
S2 1.386 + 0.027 
Q2 0.526 + 0.039 
a (rad) 2.292 + 0.013 
53 1.601 + 0.043 
q3 0.241 + 0.027 
y (rad) 1.383 + 0.010 
p (10°) 2.50 + 0.21 
Notes. HJD’ = HJD — 2450000. 


We note that the source radius ratio p2/pı = 1/46 and the source 
flux ratio gr = 0.15 would lead to inconsistent estimates of 6g, 
which would argue against the physical plausibility of this solu- 
tion. However, we do not pursue this issue because this solution 
will prove to be heavily disfavored by y?. See Sect. 3.3. 


3.3. 3L1S interpretation 


We additionally searched for a model under the 3L1S interpre- 
tation. Similar to the case of the 2L2S modeling, the 3L1S mod- 
eling was done based on the 2L1S solution, from which we 
adopted the initial parameters of (tọ, Uo, te, 52, 92, œ, p). We then 
conducted preliminary searches for the parameters related to M3, 
that is, (s3, q3, Y), via a grid approach, and finally refined the so- 
lution by letting all parameters vary. Figure 3 shows the Ay? map 
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Fig. 4. Lens-system configuration of the 3L1S model. The main panel 
shows the central anomaly region, and the inset displays the whole view 
of the lens system, where blue dots marked by M,, M>, and M; denote 
the positions of the individual lens components. Grey curves around 
the caustic represent the equi-magnification contours. The dotted circle 
centered at the origin of the coordinates represents the Einstein ring of 
the lens system. 


on the log s3—-log q3 parameter plane constructed from the pre- 
liminary grid search. The map shows that there exists a unique 
y minimum at (log s3, log q3) ~ (0.2, —0.6). 

It was found that the model identified under the 3L1S in- 
terpretation well described all the observed anomaly features in 
the lensing light curve. The lensing parameters of the best-fit 
3L1S solution are listed in Table 2 together with the value of 
x’ /dof. In Figure 2, we draw the model curve (solid curve drawn 
over the data points) and residual (presented in the panel la- 
beled "3L1S") from the model. The parameters related to Mz and 
M3 are (s2,q2) ~ (1.4, 0.53) and (53, q3) ~ (1.60, 0.24), respec- 
tively, indicating that the primary of the lens is accompanied by 
two lower-mass companions lying in the vicinity of the Einstein 
ring. If the separations among the lens components were intrinsic 
(3-dimensional), the three-body lens system would be dynami- 
cally unstable, and this indicates that one or both companions are 
likely to be aligned by chance with the primary. We found that 
it was difficult to detect microlens-parallax effects (Gould 1992) 
due to the relatively short time scale, tg ~ 15 days, of the event 
compared to the orbital period of Earth. 

The lens-system configuration of the 3L1S model is shown 
in Figure 4, in which the main panel shows the central anomaly 
region and the inset shows the whole view of the lens system. 
The coordinates of the configuration is centered at the effective 
position of the Mı—M3 pair defined by Di Stefano & Mao (1996) 
and An & Han (2002), and the grey curves encompassing the 
caustic represent the equi-magnification contours. One finds that 
the caustic induced by M3 is similar to that of the 2L1S solu- 
tion, but the 3L1S caustic and the resulting magnification pattern 
around the caustic are distorted by M3, which makes the central 
part of the caustic nested and self-intersected (Gaudi et al. 1998; 
Danék & Heyrovsky 2015a,b, 2019). According to the 3L1S so- 
lution, the anomaly appearing around f) was produced by the 
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Fig. 5. Cumulative distribution of the y* difference, Ay? = 5, 95—YX 119° 
between the 2L2S (orbital) and 3L1S models in the region of the major 
anomaly. The light curve in the upper panel is shown to indicate the 
epochs of major y? differences. 


source approach close to the tip of the swallow-tail caustic part 
induced by M3. 

From the comparison of the 3L1S and 2L2S (orbital) mod- 
els, we find that the anomaly features in the lensing light curve 
are better explained by the 3L1S interpretation than by the 2L2S 
interpretation. In Figure 5, we present the cumulative distribu- 
tion of the y? difference, Ay? = XS -= Xi ş» between the two 
models in the region of the anomaly for the detailed comparison 
of the two models. The major differences occur around tọ and 
tı, at around which Sz and S, first crossed the caustic accord- 
ing to the 2L2S model, respectively. Besides these regions, the 
3L1S model well explains the KMTA data point at the epoch of 
HJD’ = 9373.272, while the residual of this data point from the 
2L2S model, AJ = 0.656 mag, is very big. As a whole, the 3L1S 
model yields a better fit than the 2L2S model by Ay” = 53.8, 
despite the fact that the 3L1S model has 4 fewer parameters. 


4. Source star and Einstein radius 


We specify the source of the event not only for the full charac- 
terization of the event but also for the estimation of the angular 
Einstein radius. The Einstein radius is determined from the an- 
gular radius of the source, 6,, by the relation 


(1) 


where @, is deduced from the source type, and the normalized 
source radius p is measured from the modeling by analyzing the 
caustic-crossing parts of the light curve. 

We determined the source type by measuring the extinction- 
and reddening-corrected color and brightness, (V — Z, I)o,s. For 
this, we utilize the Yoo et al. (2004) method, in which the source 
is located in the instrumental color-magnitude diagram (CMD), 
and then its color and magnitude are calibrated using the centroid 
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Fig. 6. Position of the source (blue dot) with respect to the centroid 
of red giant clump (RGC, red dot) in the instrumental color-magnitude 
diagram of stars lying in the vicinity of the source. 


of the red giant clump (RGC) as a reference. The RGC centroid 
can be used as a reference for calibration because its de-reddened 
color and magnitude are known. 

Figure 6 shows the position of the source in the instrumen- 
tal CMD constructed from the pyDIA photometry of stars lying 
around the source. We note that the CMD shows only bright stars 
due to the severe extinction, A; = 3.3, toward the field. The mea- 
sured instrumental color and magnitude of the source are 


(V - I, Ds = (3.501 + 0.056, 21.080 + 0.037), (2) 


where the J and V-band magnitudes were estimated by regress- 
ing the KMTC light curve data of the individual passbands pro- 
cessed using the same pyDIA code with respect to the model. 
By measuring the offsets in the color and magnitude, A(V — J, J), 
of the source from the RGC centroid, with (V — I, D roc = 
(3.669, 18.025), together with the known de-reddened color and 
magnitude of the RGC centroid, (V — I, Dorae = (1.060, 14.377) 
(Bensby et al. 2013; Nataf et al. 2013), the de-reddened source 
color and magnitude were measured as 


(V -I, Dos = (V - I, Doac + AV - 1,1) 


3 
= (0.892 + 0.056, 17.432 + 0.037). G) 


The estimated color and magnitude indicate that the source is a 
subgiant star of an early K spectral type. 

The angular source radius was deduced from the mea- 
sured source color and magnitude. For this, we first converted 
V — I color into V — K color using the color-color rela- 
tion of Bessell & Brett (1988), and then estimated the angular 
source radius using the relation between (V — K, V) and 6, of 
Kervella et al. (2004). This procedure yielded the angular radii 
of the source and Einstein ring of 


0. = (1.26 + 0.11) yas. (4) 
and 

0. 
Og = = = (0.50 0.06) mas, (5) 
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respectively. Together with the measured event time scale, the 
relative lens-source proper motion was estimated as 


u= 2 = (12.49 + 1.53) mas yr`!. 
E 


(6) 


5. Physical lens parameters 


In this section, we estimate the physical lens parameters of the 
mass M and distance D; to the lens. These parameters can 
be sometimes constrained by measuring lensing observables of 
(tg, Og, me), where mg denotes the microlens parallax. The first 
two observables are related to the physical parameters by 


Og 
tg = —; 
H 


Og = (KMT), (7) 


where x = 4G/(c*AU) = 8.14 mas/Mo and aye) = AL — Ts = 
AU(1/D, — 1/Ds) represents the relative lens-source parallax. 
With the additional measurement of the observable 7g, the phys- 
ical parameters can be uniquely constrained by the relations 
(Gould 2000) 


au 
M = we DL = 
KTE 


ngebe + Ts ; (8) 
For KMT-2021-BLG-1122, the observables tg and Og were se- 
curely measured, but the other observable zg could not be mea- 
sured, and thus we estimate the physical lens parameters by con- 
ducting a Bayesian analysis using a Galactic model together with 
the constraints provided by the measured observables. 

In the first step of the Bayesian analysis, a large number of 
artificial lensing events were produced from a Monte Carlo sim- 
ulation. For each simulated event, we derived the locations of 
the lens and source and their transverse velocity from a Galactic 
model, and assigned the lens mass from a mass-function model. 
In the Monte Carlo simulation, we adopted the Galactic model of 
Jung et al. (2021), whose detailed description are given therein. 
With the simulated events, we then constructed the posteriors of 
the lens mass and distance by giving a weight to each simulated 
event of 


(tei — te)’ (Oei — Op)? 
i = exp(—x;/2); oe r 
w; = exp(—y; /2) Xi ot] [o(6£)] 


where (tg, @g) and [o(tge), o(0g)] represent the measured val- 
ues of the observables and their uncertainties, respectively, and 
(tE,i, 9¢,;) indicate the observables of each simulated event. 

Figure 7 shows the posteriors of the primary lens mass (top 
panel), and distances to the lens (middle panel) and source (bot- 
tom panel). We additionally present the Ds posterior to show 
the relative location of the lens and source. For each posterior, 
we mark the contributions by the disk and bulge lens popula- 
tions with curves drawn in blue and red, respectively. The masses 
of the individual lens components estimated from the Bayesian 
analysis are 


(9) 


Mı = 0.477033 Mo, (10) 
My = 0.24*313 Mo, (11) 
and 

Mz = 0.117906 Mo. (12) 
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Fig. 7. Posteriors of the primary lens mass (top panel) and distance to 
the lens (middle panel) and source (bottom panel). In each panel, the 
blue and red curves represent the contributions by the disk and bulge 
lens populations, respectively, and the black curve is sum of the two 
lens populations. The solid vertical line indicates the median value, and 
the two dotted lines represent the 10 range of the distribution. 


The estimated masses of all the lens components are in the mass 
regime of M dwarfs, and thus the lens is a triple system com- 
posed of three low-mass stars. The estimated distance to the lens 
system is 

Di = 657-55, kpc, (13) 
and the probabilities for the lens to be in the disk and bulge are 
33% and 67%, respecively. For each parameter, we present the 
median as a representative value and the lower and upper limits 
are estimated as the 16% and 84% ranges of the posterior distri- 
bution. The projected separations of the secondary and tertiary 
lens components from the primary are 


412 = s2D10g = 3.477) 1) AU, (14) 
and 
13 = 53D, 0g = 4.014383 AU, (15) 


respectively. As noted, the fact that the projected separations a, 2 
and a, 3 are similar to each other suggests that one or both lens 
companions are aligned by chance with the primary. 

We note that KMT-2021-BLG-1122L is the first triple stellar 
system detected via microlensing. In total, there exist 12 con- 
firmed 3L1S events: OGLE-2006-BLG-109 (Gaudi et al. 2008; 
Bennett et al. 2010), OGLE-2012-BLG-0026 (Han et al. 2013), 
OGLE-2018-BLG-1011 (Han et al. 2019), OGLE-2019-BLG- 
0468 (Han etal. 2022d), OGLE-2021-BLG-1077 (Han et al. 
2022a), OGLE-2006-BLG-284 (Bennett et al. 2020), OGLE- 
2007-BLG-349 (Bennett et al. 2016), OGLE-2008-BLG-092 
(Poleski et al. 2014), OGLE-2013-BLG-0341 (Gould et al. 
2014), OGLE-2018-BLG-1700 (Han et al. 2020b), KMT-2019- 
BLG-1715 (Han etal. 2021a), and KMT-2020-BLG-0414 
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(Zang et al. 2021).* Among them, the lenses of the first 5 events 
are two-planet systems, and the lenses of the other 7 events are 
binary-star systems possessing planets. Therefore, the common 
trait of the previously discovered microlensing triple systems is 
that at least one of the lens components is a planet. Two planets 
with similar separations from a host can be dynamically stable 
if the planets are in mean-motion resonance (Madsen & Zhu 
2019). A planet in a binary system can be dynamically stable 
if the planet orbits one star of a wide binary stellar system or 
the barycenter of a closely-spaced binary system. Microlens- 
ing detection of a triple stellar system is difficult because a 
system with similar masses and intrinsic separations among 
the components would be dynamically unstable. One way 
for such a system to be detected via microlensing is that the 
companions of the triple system are closely aligned with the 
primary so that they lie around the Einstein ring of the primary. 
This alignment requires projection effect of a low probability, 
and this explains the relative rareness of triple stellar systems 
detected by microlensing. 


6. Summary and conclusion 


We presented the analysis of the microlensing event KMT-2021- 
BLG-1122, which was investigated as a part of the project to 
analyze anomalous lensing events in the previous KMTNet data 
with no suggested plausible models. We confirmed that the light 
curve, characterized by three major anomaly features, of the 
event could not be explained with the usual 2L1S or 1L2S in- 
terpretations, but it was found that a 2L1S solution obtained 
from the modeling with the exclusion of the data around the sec- 
ond anomaly feature could approximately describe the other two 
caustic-crossing features. 


We found that all the anomaly features could be explained 
by introducing a tertiary lens component, and thus the lens 
is a triple system. The masses of the individual lens compo- 
nents estimated from the Bayesian analysis are (M1, Mz, M3) ~ 
(0.47 Mo, 0.24 Mo, 0.11 Mo), and the companions are separated 
in projection from the primary by (a1,2,a1,3) ~ (3.5,4.0) AU. 
The estimated masses of all the lens components are in the mass 
regime of M dwarfs, and thus the lens is a triple system com- 
posed of three low-mass stars. This is the first triple stellar-lens 
system detected via microlensing. 
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+ Besides these confirmed cases, there are 6 candidate triple-lensing 
events, including OGLE-2014-BLG-1722 (Suzuki et al. 2018), OGLE- 
2018-BLG-0532 (Ryu et al. 2020), KMT-2019-BLG-1953 (Han et al. 
2022a), KMT-2019-BLG-0304 (Han et al. 2021b), OGLE-2019-BLG- 
1470 (Kuang etal. 2022), and KMT-2021-BLG-0240 (Han et al. 
2022b). For these events, 3L1S interpretations yielded best-fit models, 
but the signals of the tertiary lens components were not firmly con- 
firmed either due to the subtlety of the signals or the degeneracy with 
other interpretations. Under the 3L1S interpretations of these events, the 
lenses belong to either the category of two-planet systems or the cate- 
gory of planets in binaries similar to the cases of the confirmed 3L1S 
events. 
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